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摘 要 : 地 表 履 被 作为 自然 过 程 和 人 类 活动 共同 作用 的 重要 地 表 景 观 特征 ,对 全 球 或 局 地 气候 、 水 


热 循环 、 物 质 传输 及 陆 


4 生态 系统 多 样 性 等 影响 深远 。 利 用 年 内 时 序 遥 感 影 像 自动 提取 不 同 地 表 
履 被 类 型 的 方法 ,以 新 疆 乍 康 地 区 为 研究 目标 ,组 织 2016 年 植被 全 生长 季 的 Landsat 8 OLI 地 表 反 


射 时 序 影 像 ,研究 不 同 物候 期 植被 冠 层 的 纹理 响应 信息 ,考察 研究 区 典型 地 表 履 被 类 型 在 3 一 11 月 
多 波段 波谱 、 归 一 化 植被 指数 (NDVI) 及 增强 型 植被 指数 (EVI) 的 时 序 特征 ,构建 提取 地 表 覆 被 类 
型 的 策略 规则 ,形成 时 序 特征 匹配 方法 ,将 其 应 用 于 2018 年 研究 区 地 表 履 被 填 图 的 提取 。 最 后 , 基 


Æ 
于 高 分 辨 率 卫 星 影 像 和 野外 实地 调查 对 随机 选取 的 2500 个 样 点 进行 对 比 验证 。 结 果 表 明 : 提 取 结 
果 和 验证 数据 一 致 性 较 好 ,总 体 精 度 为 97.29% ,Kappa 系 数 为 0.9655, 且 实地 考查 结果 显示 本 方法 在 


复 播 作 物 识别 和 有 效 降 低 单 一 时 相 中 "异物 同 谱 " 现 象 上 展示 潜在 优势 。 
关键 词 : 时 序 遥 感 数据 ; 植被 物候 ; 光谱 特征 ; 归 一 化 植被 指数 ; 增强 型 植被 指数 ; 地 表 覆 被 填 图 
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地 表 有 覆 被 逐渐 成 为 气候 变化 .国情 监测 .生态 
评价 和 可 持续 发 展 规划 研究 中 的 重要 变量 。 近 年 
来 ,虽然 已 有 不 少 关于 全 球 地 表 履 被 数据 用 于 满足 
生物 多 样 性 和 农业 管理 的 研究 ,但 是 数据 适用 性 仍 
显 不 足 悦 。 不 仅 全 球 约 60% 的 农田 制图 精度 仍 需 提 
高 中 ,粮食 估 产 和 农 情 监测 对 此 也 有 更 高 精度 的 需 
Ko 虽然 地 表 覆 被 产品 分 类 精度 逐步 提升 ,但 地 表 
窗 被 作为 地 球 表面 物质 和 自然 属性 的 综合 体 中 , 受 
气候 环境 等 影响 ,尤其 是 在 地 表 蔡 散 严 重 的 干旱 
区 ,表现 出 很 强 的 脆弱 性 外 ,其 识别 精度 尚未 满足 广 
大 用 户 日 益 增 长 的 应 用 需求 。 

随 着 地 表 分 类 精度 需求 的 提高 和 数据 的 累积 ， 
对 地 数据 逐渐 从 零星 到 系统 的 信号 采集 和 组 织 模 
式 。 早 期 利用 单一 时 相 遥 感 影像 进行 地 表 履 被 类 
型 的 识别 方法 简便 快捷 ,但 就 地 物 信息 提取 仅 做 时 
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间 维 上 有 限 离散 甚至 孤立 的 采样 ,难以 捕捉 许多 自 
然 地 表 随 时 间 推 移 或 者 季节 更 蔡 的 渐进 变化 特点 ， 
潜在 缩减 了 对 “ 同 谱 异 物 ” 正 确 解 译 的 机 会 ,常常 不 
易 满 足 动态 监测 精度 要 求 。 而 加 密 的 时 间 序 列 遥 
感 影像 则 可 以 有 效 采 集 地 表 覆 被 随 物 候 或 季 相 变 
化 的 节律 特点 ,为 进一步 细 化 地 表 类 别 和 提高 不 
同 地 表 覆 被 识别 精度 提供 可 能 。 

综合 地 物 时 序 光 谱 特 征 或 植被 指数 信息 提取 
地 表 履 被 正在 逐渐 被 越 来 越 多 的 学 者 关注 。 如 利 
用 多 时 相 Landsat 8 OLI 数 据 深度 挖掘 地 物 的 光谱 特 
征 及 归 一 化 植被 指数 (Normalized difference vegeta- 
tion index,NDVI) 的 时 序 信息 ,选取 某 些 关键 时 相 构建 
不 同 作物 的 决策 树 模 型 ,提取 作物 种 植 结构 信息 ””。 
胡 勇 等 "根据 训练 样本 和 参考 影像 光谱 特征 将 其 
在 时 间 维 进行 扩展 以 提高 分 类 精度 。 王 新 洁 等 
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基于 Landsat 等 中 分 状 率 遥感 数据 定量 提取 地 物 光 
谱 和 时 序 特 征 ,同时 应 用 各 种 分 类 器 及 多 源 算法 ， 
将 地 物 分 为 裸 地 ,耕地 、 林 地 等 基础 类 别 ,生产 连续 
地 表 覆 被 填 图 。 但 植被 物候 信息 在 中 高 分 辨 率 地 
表 覆 被 提取 中 的 应 用 仍 很 缺乏 ,尤其 在 农 情 管理 、 
智能 农业 或 精准 农业 监测 等 领域 沿 需 更 为 精细 的 
农作物 类 别 填 图 。 因 此 ,如 何 利 用 时 间 序 列 遥 感 影 
像 定量 描述 农作物 的 物候 变化 规律 ,表征 不 同 地 表 
窗 被 的 特点 ,提高 地 表 窗 被 精度 具有 重要 意义 。 
综 上 所 述 ,时 间 序 列 方法 在 多 种 地 物 分 类 尤其 
是 农作物 类 别提 取 中 有 待 进一步 加 强 。 本 文 基于 
作物 生长 季 时 序 的 中 分 辨 率 地 表 反 射 数 据 ,提取 典 
型 地 物 多 波段 波谱 和 NDVI 及 增强 型 植被 指数 (En- 
hanced vegetation index , EVI) 的 时 序 特 征 ,建立 特征 
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规则 提取 实验 区 作物 类 别 信息 。 
材料 与 方法 


1.1 研究 区 概况 

本 文选 取 西 部 绿洲 城市 3 旱 康 地 区 (87?53'~ 
88°4'E ,44°15'~44°22'N, R D AUK ,其 位 于 新 
疆 维吾尔 自治 区 昌吉 回族 自治 州 ,地 处 天 山东 段 北 
FE 、 准 噶 尔 盆地 的 南 缘 , 中 部 为 由 天 山 融 雪 河 流 冲 
积 而 成 的 平原 ,是 绿洲 农业 的 集中 分 布 区 ,具有 多 
样 的 作物 类 型 。 新 疆 旱 康 地 区 属于 典型 的 温带 大 
陆 性 干旱 气候 ,水 资源 荐 乏 , 莱 发 强烈 ,主要 灌溉 水 
源 是 天 山 冰雪 融 水 与 细小 河流 ,对 精准 农业 的 需求 
更 为 迫切 ,因此 对 该 区 域 地 表 获 被 类 型 的 提取 ,万 
其 是 农作物 信息 的 准确 提取 变 得 尤为 重要 。 
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图 1 研究 区 真 彩色 影像 及 地 理 位 置 


Fig. 1 True color images and geographical location of the study area 


1.2 研究 方法 

1.2.1 时 序 特 征 匹 权 不 同 地 表 覆 被 (尤其 是 农 作 
物 ) 生 育 周期 内 的 地 表 反 射 波谱 及 植被 指数 各 具 特 
色 ,特征 明显 。 如 图 2 为 不 同类 型 植被 生长 发 育 特 
征 随时 序 典 型 分 异 情况 。 一 些 秋末 初冬 播种 的 作 
物 , 如 小 麦 , 在 春季 进入 发 育 旺 盛 时 期 ,到 夏 初 就 发 
育成 熟 进入 收割 期 ; 而 春播 作物 出 苗 就 已 5.6 月 份 ， 
经 过 各 具 特 色 的 发 育 过程 , 到 夏 未 或 秋季 结束 生长 
进入 收割 期 ,如 玉米 .棉花 等 作物 ;此 外 城市 绿化 、 
乡村 绿 植 ,防护林 或 防风 带 等 也 各 有 其 独特 的 生长 
发 育 和 休 眼 特征 。 结 合 植被 物候 期 及 纹理 特征 选 


1.01 一 -小麦 
0.8 | -棉花 
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时 间 ( 月 -日 ) 


图 2 典型 植被 的 时 间 序列 特征 图 
Fig.2 Illustration of temporal features 


implied in typical vegetation 
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FHEWY 


取 典 型 地 表 履 被 样本 点 ,提取 不 同 植被 在 波谱 空间 
全 生育 期 内 绿 起 成 熟 / 丰 叶 、 调 零 /收割 多 波段 波谱 
和 植被 因子 的 渐进 变化 特征 ,组 成 特征 向 量 
x, Expt) ,其 中 vi 表示 第 i 种 植被 ; x(j= 
1,2,…,n) 表示 所 提取 的 n 维 时 序 特征 ,用 来 定量 
征 植被 v_i 在 不 同 生育 期 的 典型 特征 值 。 

以 欧 氏 距离 中 和 最 近邻 原则 为 主 判 别 待 识 别 
像 元 Vp 的 精细 类 别 归属 。 首 先 依据 公式 (1) 估 算 
Vp 的 时 序 特征 向 量 (yyy) 与 m 个 先 验 特征 样 
本 (xpa) 的 特征 距离 (D, 1,D,,,…,D, n) ,再 根 
据 最 近邻 原则 进行 类 别 匹配 ,获得 最 大 似 然 的 精细 
类 别 。 


D,= (1) 


È =y) 
UP: D, 指 特征 距离 ,其 中 k=1,2,…,m,jell,+%); 
xn 指 训练 样本 时 间 序 列 的 特征 值 ; y, 指 待 分 类 样本 
的 时 间 序 列 特征 值 。 

1.2.2 技术 流程 图 3 所 示 为 主要 技术 流程 。 首 先 
Xf Landsat 8 OLI 数 据 进行 预 处 理 (模块 A) ,为 有 效 
提取 地 表 履 被 纹理 .波谱 特征 服务 ;融合 物候 信息 ， 
构建 典型 地 表 覆 被 的 先 验 特征 向 量 (模块 B) ,为 自 
动 识别 建立 规则 ;再 将 其 应 用 到 待 识别 影像 序列 ， 
自动 进行 特征 匹配 ,获取 最 似 然 的 精细 地 物 类 别 
(模块 C); 最 后 与 验证 数据 进行 对 比 与 精度 分 析 ( 模 
块 D)。 

1.3 数据 来 源 

1.3.1 物候 数据 新 疆 旱 康 地 区 地 表 观 测 农作物 生 
长 发 育 的 物候 期 ( 旬 / 月 ) 见 表 1。 该 地 区 典型 地 表 覆 
被 在 生长 发 育 过 程 中 ,会 呈现 不 同 的 特征 。 如 小 麦 


o OOO‘ 
A 图 像 预 处 理 Landsat 地 表 反 射 率 产品 
BERRIE yy. 
sisal 图 像 融合 
DX Val 
十 比 验证 图 像 裁剪 


| 典型 植被 物候 信息 

典型 地 物 纹理 特征 分 析 

: 典型 地 物 多 

人 波段 反射 率 提取 

| 典型 多 光谱 

| 时 序数 据 植被 指数 时 序 

i's 特征 分 析 与 描述 

有 

IE a 

| eae 时 序 特征 匹配 | 

h 相 和 遥感 数据 计算 最 小 欧 氏 距离 | 

人 | 
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ooo 间接 对 比 | 

:| 对比 验证 数据 集 定性 验证 

D e 
vy 


综合 分 析 与 讨论 | 


图 3 信息 处 理 流 程 


Fig. 3 Flow chart of information process 


播种 早期 不 覆 膜 ,9 月 中 下 旬 播 种 ,翌年 6 月 中 下 旬 
成 熟 ;而 棉花 和 玉米 播种 早期 均 覆 膜 , 棉 花 每 年 4 月 
播种 ,11 月 中 旬 收 割 结束 ,玉米 每 年 4 月 播种 ,9 月 
底 收 割 ,其 收获 期 间 绿叶 仍 可 能 存在 ;水 稻 在 分 蓝 拔 
节 期 需 灌 入 大 量 的 水 ; 油 葵 生育 期 较 短 ,持续 120 d 
左右 ,其 在 4 月 末 播 种 。 


表 1 新 疆 阜 康 地 区 农作物 物候 期 
Tab.1 Phenological period of crops in Fukang City, Xinjiang 


农作物 生育 期 RG Sh BE DRAKE 返青 一 拔节 拔节 一 抽穗 抽穗 一 成 熟 
小 麦 下 /9 一 上 /11 上 /1 一 下 /3 下 /3 一 中 /4 中 /4 一 中 /5 中 /5 一 下 /6 
农作物 生育 期 播种 一 出 苗 出 苗 一 拔节 拔节 一 抽 惩 抽穗 一 乳 熟 乳 熟 一 成 熟 
玉米 下 /4 一 中 /5 中 /5 一 下 /6 下 /6 一 下 /7 下 /7 一 中 /8 中 /8 一 中 /9 
农作物 生育 期 Gi h th Hi BF MIF AE FREER TEUS 
棉花 下 /4 一 中 /5 中 /5 一 中 /6 中 /6 一 下 /7 /7 一 下 /9 下 /9 一 上 /11 
农作物 生育 期 RBS ee 拔节 一 抽穗 HEEE 灌浆 一 成 熟 
水 稻 下 /4 一 中 /5 中 /5 一 下 /6 下 /6 一 上 /8 上 /8 一 上 /9 上 /9 一 下 /10 
农作物 生育 期 Pi i tH iF IEE 开花 一 成 熟 

Pips 下 /4 一 中 /5 中 /5 一 上 /6 上 /6 一 下 /7 下 /7 一 下 /8 


注 : 下 /9 等 表示 农作物 生长 发 育 的 物候 期 ( 旬 / 月 )。 
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1.3.2 验证 数据 ”本文 选取 覆盖 研究 区 的 高 分 辨 京 
开放 遥感 影像 (谷歌 地 图 Google Earth ) 和 野外 实地 
调查 数据 作为 验证 数据 。 随 机 选取 2500 个 样本 点 ， 


小 麦 (图 4a),7.8 月 的 油 芋 (图 4b) 和 玉米 (图 4c) 等 ; 
而 近 红 外 波段 B5(0.85~0.88 um ) 随 植被 出 苗 拔节 等 
生长 变化 其 冠 层 绿色 叶片 增加 , 近 红 外 波段 反射 率 


包括 小 麦 .玉米 SS 水稻、 棉花 、 葡 萄 、 居 民 地 E 
地 、 裸 土 和 复 播 10 类 地 物 ,利用 谷歌 地 图 对 其 进行 
验证 ,同时 以 2019 年 5 月 18 日 野外 实地 调查 的 葡 
萄 .水稻 .小 麦 和 油 蕉 作为 辅助 数据 来 进行 验证 。 
13.3 影像 数据 及 预 处 理 “本文 所 采用 的 遥感 影像 
数据 来 源 于 美国 地 质 调 查 局 (https://earthexplorer. 
usgs.gov/) 所 发 布 的 已 做 过 地 理 定 标 ,大气 校正 的 地 
表 反 射 率 数 据 产 品 汪 和 ,可 有 效 反 应 地 表 波 谱 响 应 
音 息 。 经 对 比 筛 选 ,选取 2016 年 (14 个 时 相 ) 和 
2018 年 (12 个 时 相 ) 云 量 较 少 ,植被 生长 季 时 相 较 多 
的 Landsat 8 OLI 数 据 。 其 中 ,2016 年 的 14 个 时 相 作 
为 提取 波谱 及 植被 指数 特征 的 数据 ,2018 年 的 12 个 
时 相 数 据 作为 地 表 禾 被 分 类 数据 。 

为 满足 实验 需求 ,对 研究 区 遥感 数据 进行 了 预 
处 理 操作 。 将 遥感 影像 进行 4.3.2 波 段 的 融合 形成 
真 彩色 图 像 ,并 对 息 康 地 区 农业 集中 分 布 研究 区 进 
行 裁剪 操作 ,图 la 为 旱 康 市 真 彩 色 图 像 , 图 lb 为 所 
裁剪 的 研究 区 真 彩色 图 像 。 再 进行 研究 区 时 相 数 
据 间 的 配 准 校 验 ,保证 不 同时 期 获取 的 影像 精确 配 
准 。 最 后 ,对 地 面 观测 农作物 物候 期 与 遥感 影像 数 
据 采 集 日 期 归 一 化 整合 处 理 。 此 外 ,准备 调查 
Google Ear 由 在 该 地 区 清晰 影像 和 整理 实地 调查 记 
录 ,为 分 类 验证 做 准备 。 


2 结果 与 分 析 


2.1 时 间 序 列 多 波段 波谱 特征 分 析 

不 同 地 表 覆 被 的 反射 波谱 在 不 同 季 相 ,尤其 是 
植被 ,天然 拥 有 各 具 特 色 的 生长 发 育 周期 ,大 都 差 
异 鲜明 。 根 据 旱 康 农作物 产业 结构 年 鉴 数 据 和 新 
性 旱 康 地 区 农作物 物候 期 可 知 研 究 区 的 农作物 主 
要 有 小 麦 TSS .玉米 .水稻 .葡萄 和 机 花 等 。 图 4 为 
旺 康 地 区 8 种 典型 覆 被 类 型 样本 点 2016 年 Landsat 
8 0LI 中 16 个 时 相 多 光谱 时 序 变 化 曲线 。 总 体 上 ， 
可 见 光 波段 包括 蓝 波段 B2(0.45~0.51 um) , 绿 波段 
B3(0.53~0.59 um) 和 红 波 段 B4(0.64~0.67 p.m) 在 未 
有 植被 履 盖 时 ,呈现 裸 土 特性 (图 币 )。 一 旦 作物 出 
苗 或 其 他 植被 返青 ,可见光 波段 会 明显 被 抑制 , 且 
绿 波段 反射 率 比 其 他 可 见 光 的 反射 率 高 ,如 4、5 月 


变化 明显 , 且 在 植被 生长 期 , 相 比 其 他 波段 反射 率 
值 较 高 ,如 6、7 月 的 玉米 (图 4c) 及 8、9 月 的 棉花 (图 
4f) 等 ;短波 红外 1B6(1.57~1.65 jm) 对 水 具有 强 吸 
收 性 , 故 在 植被 生长 期 其 反射 率 较 低 , 非 植被 生长 
期 反射 率 较 高 ,说明 土 壤 含 水 较 少 ,地 表 较 为 干 早 
(图 4h)。 且 比 其 他 地 物 ,水 稻 受 灌水 影响 ,其 各 个 
波段 的 反射 率 较 低 (图 4d)。 
2.2 时 间 序 列 NDVI 特 征 分 析 

研究 表明 NDVI 具 有 捕获 植被 生长 状况 的 能 
力 。 不 同 植被 的 内 在 组 织 结构 .色素 含量 和 和 冠 层 组 
分 比例 在 生长 发 育 过 程 中 不 断 变 化 , 红 光 和 近 红 外 
波段 对 此 响应 相对 敏感 。NDVI 定 义 中 [ 见 公式 (2)] 
近 红 外 和 红 光 2 个 波段 差 值 由 在 放大 植被 响应 信 
号 ,并 通过 归 一 化 增强 不 同情 况 下 取 值 的 定量 可 比 
性 ,是 探讨 植被 特征 的 有 效 因子 ""。 

NDVI = Am Puro (2) 
Pu + Prev 

式 中 : Prio 和 pw 分 别 指 红 波段 和 近 红 外 波段 的 地 
表 反 射 率 。 时 间 序 列 的 NDVI 可 以 在 同一 作物 生育 
期 内 进行 多 次 采样 ,潜在 勾勒 出 不 同 作物 生长 发 育 
进展 ,丰富 分 类 信息 ,降低 “异物 同 谱 ” 几 率 ,以 提高 
准确 提取 精度 。 一 般 情况 下 , 随 着 植被 的 快速 生 
长 , 叶 面 积 持续 增加 ,卫星 观测 的 视 场 内 信号 增强 ， 
而 红 光 和 近 红 外 波段 在 光谱 空间 对 绿色 生长 变化 
的 联动 反映 ,使 不 同 地 表 覆 被 展现 出 不 同 的 波谱 特 
征 , 综 合 分 析 图 5a 可 以 得 到 : 

(1) 小 麦 4 月 中 旬 拔 节 期 前 后 , 干 物质 积累 快 生 
长 旺盛 ,NDVI 曲 线 迅 速 上 升 直至 5 月 中 下 名 抽 穗 期 ， 
NDVI 值 达到 峰值 ( 约 0.8) ,保持 约 2 周 的 饱和 状态 ， 
随 着 小 麦 灌 浆 成 熟 逐 渐 降 低 , 进 入 7 月 中 旬 收 制 期 。 
若 复 播 后 期 呈现 种 植 作 物 生 育 期 的 特征 ( 见 复 播 曲 
线 ) , 若 无 复 播 NDVI 向 裸 土 特征 逼近 , 取 值 约 为 0.1。 

(2) 油 蕉 7 月 开花 ,NDVI 值 达到 峰值 ( 约 0.7) ,8 
月 下 旬 油 给 成 熟 收割 ,NDVI 值 基本 保持 在 0.1 左 右 。 

(3) 玉米 和 油 葵 的 生育 规律 较为 相似 ,但 玉米 6 
月 下 旬 处 于 拔节 期 ,其 NDVI 值 达到 峰值 ( 约 0.8)， 
且 玉 米 7 月 中 旬 到 8 月 上 旬 处 于 抽雄 期 ,玉米 的 雄 
花 呈 淡 黄 色 穗 状 ,NDVI 值 会 呈现 低谷 ,NDVI 值 由 
0.87 下 降 至 0.83 ,又 上 升 到 0.86。 
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主 :B1 表示 海岸 波段 ;B2 表 示 蓝 波段 ;B3 表示 绿 波 段 ;B4 表 示 红 波段 ;B5 表 示 近 红外 波段 ;B6 表 示 短 波 红外 1;B7 表 示 短 波 红外 2。 


图 4 旱 康 典型 地 物 的 多 光谱 时 序 变化 


Fig.4 Dynamic changes of multispectral over typical land covers of Fukang City 
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注 :NDVIEVI 分 别 表示 归 一 化 植被 指数 .增强 型 植被 指数 。 


图 5 阜 康 典 型 地 物 的 NDVIEVI 时 序 变 化 
Fig. 5 Dynamic changes of NDVI and EVI over typical land covers of Fukang City 


(4) 水 稻 5 月 中 旬 ab FABER TTH , 需 向 稻田 
中 灌 入 大 量 的 水 ,此 时 NDVI 值 达到 低谷 ,4 月 上 多 
NDVI 值 为 0.18 ,灌水 后 ,5 月 中 旬 NDVI 值 下 降 至 
0.08,5 月 下 旬 NDVI 值 又 上 升 到 0.18。 

(5) 棉花 9 月 处 于 吐 祭 期 ,其 NDVI 值 基本 维持 
为 0.8 左 右 , 而 此 时 其 他 作物 逐渐 收割 ,NDVI 曲线 
呈现 下降 趋势 。 

(6) 葡萄 是 多 年 生 果树 ,4 月 温度 上 升 后 生长 发 
育 ,到 6 月 下 旬 NDVI 达 到 峰值 ( 约 0.7) ,随时 间 推 移 
开始 结果 实 , 由 于 葡萄 果实 呈现 紫色 ,NDVI 下 降 至 
0.6 左 右 , 大 约 需 35~50 d 果 实 成 熟 ,期 间 其 他 幼 果 也 
逐渐 成 熟 ,持续 到 10 月 , 待 果 实 摘 掉 后 , 仍 有 绿叶 存 
在 NDVI 曲 线 呈 现 上 升 的 趋势 ,10 月 到 11 月 进入 休 
眠 期 其 呈 下 降 趋 势 。 

(7) 裸 土 的 NDVI 值 一 般 较 低 , 基 本 保持 在 0.1 
左右 ,其 NDVI 值 一 般 不 随 着 时 间 推 移 而 改变 。 

(8) 草地 在 温度 适宜 时 ,进行 生长 ,5 月 中 旬 到 
10 月 上 旬 , 其 NDVI 值 维持 在 0.2~0.3 之 间 。 

(9) 居民 地 包括 建筑 物 ,以 及 与 居住 相关 的 
其 他 生活 设施 (如 道路 .园林 绿化 等 ) ,经 过 考察 得 
知 建筑 设施 周围 存在 树木 、 人 工 绿化 等 植被 ,一 般 
居民 地 为 主 ,植被 较 少 ,其 NDVI 值 为 0.2 左 右 ,与 
草地 的 时 序 NDVI 曲线 较为 相似 。 
2.3 ”时间 序列 EVI 特 征 分 析 

NDVI 高 端 取 值 时 存在 饱和 现象 ,考虑 EVI 
[ 见 公式 (3)] 辅 助 分 析 ,在 一 定 程 度 上 克服 NDVI 高 
端 取 值 的 饱和 现象 , 增 大 稠密 绿叶 植被 间 的 差异 
性 ,提高 丰 叶 期 不 同 植被 识别 的 敏感 性 1。 

2.5 X (Pum Pren) (3) 
Print Orn — 7-SP giue t 1 


Ne 


EVI = 


式 中 : pw 指 蓝 波段 的 地 表 反 射 率 。 如 图 5b 所 示 
为 研究 区 典型 地 物 时 序 EVI 曲 线 , 可 以 得 出 :5 月 中 
旬 到 5 月 末 , 小 麦 的 NDVI 值 变化 较 小 而 EVI 值 在 增 
加 ,玉米 在 抽雄 期 NDVI 值 和 EVI 值 从 7 月 中 旬 至 7 
月 下 旬 逐 渐 下 降 ,7 月 下 旬 到 8 月 上 旬 逐 渐 上 升 ,但 
EVI 值 的 变化 幅度 较 大 ;水 稻 在 分 碍 拔节 期 需 向 稳 
田中 灌 和 人 大 量 的 水 , 故 其 NDVI 值 和 EVI 值 4 月 上 名 
到 5 月 中 旬 逐 渐 下 降 , 而 5 月 中 旬 到 5 月 下 旬 逐 渐 上 
升 , 且 EVI 值 的 变化 幅度 较 大 ; 油 葵 .葡萄 .棉花 的 
EVI 值 每 个 时 间 节 点 均 呈 现 一 定 变 化 的 趋势 ;草地 
和 居民 地 的 时 序 NDVI 曲 线 特 征 较为 相似 , 而 时 序 
EVI 曲 线 较为 不 同 ,可 说 明 2 种 地 物 的 差异 性 ; 裸 土 
的 时 序 EVI 曲 线 变化 较为 平缓。 


3 讨论 
3.1 方法 应 用 


根据 上 述 基 于 时 序 光 谱 及 植被 指数 特征 构建 
的 方法 自动 识别 地 表 覆 被 情况 。 综 合 分 析 时 间 序 
列 NDVI 曲 线 和 时 间 序 列 EVI 曲 线 的 变化 可 以 得 知 ， 
NDVI 有 曲线 和 EVI 曲 线 趋势 相同 ,其 中 NDVI 值 经 过 
归 一 化 处 理 , 具 有 定量 可 比 性 ,选用 时 间 序 列 NDVI 
值 对 地 表 覆 被 信息 进行 提取 ,EVI 值 作为 辅助 信 
息 。 本 文 将 已 提取 的 2016 年 时 间 序 列 地 表 覆 被 的 
NDVI 值 作为 匹配 模板 ,将 2018 年 时 间 序 列 地 表 覆 
被 的 NDVI 值 作为 待 分 类 样本 。 基 于 典型 样本 建立 
的 匹配 模板 与 2018 年 12 个 有 效 时 相 进 行 特 征 匹 
配 , 根 据 欧 氏 距离 和 最 近 令 原则 判定 2018 年 NDVI 
影像 的 每 个 像 元 对 应 的 类 别 ,完成 2018 年 时 康 研究 
区 的 地 表 精 细 类 别 填 图 (图 6c)。 如 图 6c 所 示 可 知 ， 
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图 6 2018 年 蜂 康 研究 区 的 地 表 覆 被 分 类 


Fig. 6 Land covers classification of Fukang in 2018 


该 研究 区 的 主要 农作物 主要 包括 小 麦 THE EK 
水 稻 葡萄 和 棉花 ,其 中 玉米 种 植 面 积 较 多 ,零星 分 
布 棉花 和 水 稻 的 种 植 面积 均 较 小 ,和 芽 稀 主要 集中 分 
布 在 西北 部 地 区 。 东 北部 地 区 的 圆 形 地 块 是 由 于 
干旱 区 农业 灌溉 技术 管理 模式 的 田 块 。 此 外 ,草地 
和 裸 土 相间 分 布 ,居民 地 包括 建筑 物 、 道 路 .园林 绿 
化 等 。 
3.2 地 表 覆 被 验证 

因 本 区 域 精细 到 作物 种 类 的 高 分 辨 率 地 表 覆 
被 填 图 未 有 开放 产品 可 做 横向 对 比 验 证 ,我 们 采取 
空间 采样 的 考察 验证 。 以 谷歌 地 图 和 新 疆 生 态 与 
地 理 研 究 所 提供 的 野外 实地 考察 照片 对 2018 年 地 
表 履 被 分 类 结果 进行 验证 。 在 2018 年 地 表 履 被 十 


图 上 随机 选取 100 个 150 mx150 m 的 验证 样本 区 
域 , 参 考 谷歌 地 图 地 表 履 被 数据 和 野外 实地 考察 昭 
片 对 其 进行 人 工 目 视 验 证 ,本文 实 地 调查 验证 的 作 
物 包 括 水 稻 和 葡萄 (图 6a .图 6b) ,由 于 这 2 种 多 年 生 
作物 在 2018 年 和 2019 年 种 植 区 域 基本 保持 不 变 ， 
而 其 他 农作物 可 能 轮作 而 使 种 植 区 域 变化 ,因此 可 
以 利用 包含 经 纬度 信息 的 实地 考查 作物 信息 进行 
分 类 精度 验证 。 综 合 以 上 验证 信息 ,采用 混 消 矩阵 
对 其 进行 精度 评价 ,根据 选 定 的 验证 样本 统计 每 种 
地 物 的 像 元 个 数 ,计算 混淆 和 矩阵 ,其 分 类 总 体 精度 
达到 97.2% ,Kappa 系 数 达 0.9655 ,计算 得 到 各 类 别 
的 生产 精度 和 用 户 精 度 见 表 2。 综 合 分 析 表 2 可 以 
得 到 ;本文 所 构建 的 提取 地 表 履 被 类 别 的 规则 在 一 
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表 2 2018 年 阜 康 研究 区 的 地 表 履 被 填 图 混淆 矩阵 


Tab. 2 Confusion matrix of surface cover mapping in Fukang research area in 2018 


地 表 真实 值 
小 麦 油茶 玉米 水 稻 葡萄 棉花， RE we 居民 地 复 播 总 数 用户 精 度 /% 
分 类 结果 

小 麦 170 0 0 0 0 0 0 0 0 0 170 100.00 
ub 0 298 7 0 1 0 0 0 0 0 306 97.39 
E 0 9 660 4 4 1 0 0 0 0 678 97.35 
水 稻 0 0 4 39 0 0 0 0 0 0 43 90.70 
葡萄 0 0 3 1 175 2 0 0 0 0 181 96.69 
棉花 0 0 2 0 2 45 0 0 0 0 49 91.84 

生 0 0 0 0 0 0 265 7 3 0 275 96.36 
草地 0 0 0 0 0 0 7 677 2 0 686 98.69 
居民 地 0 0 0 0 0 0 5 5 93 0 103 90.29 
复 播 1 0 0 0 0 0 0 0 0 8 9 88.89 
总 数 171 307 676 44 182 48 277 689 98 8 总 体 精度 =97.2% 


生产 者 精度 /% 99.41 97.07 97.63 88.64 96.15 


定 程度 上 降低 "异物 同 谱 ” 概率, 提高 了 对 地 物 的 识 
别 精度 。 其 中 小 麦 的 时 间 序 列 曲 线 特征 明显 ,其 物 
候 期 与 其 他 作物 差异 较 大 , 故 其 分 类 精度 较 高 ;水 
稳 和 玉米 虽 在 5 月 中 旬 存 在 较 大 差异 ,但 其 他 时 相 
特征 规律 较为 相似 ,葡萄 和 棉花 时 序 特征 变化 趋势 
较为 相近 , 且 由 于 各 种 作物 长 势 的 区 别 ,每 块 农田 
受降 雨 温度 等 影响 ,使 得 在 一 定 程 度 上 造成 误 分 ; 
裸 土 .草地 和 居民 地 三 者 曲线 变化 趋势 较为 相似 ， 
故 在 一 定 程 度 上 造成 了 误 分 ; 复 播 和 其 他 地 物 时 序 
特征 有 较 大 差异 , 故 其 提取 精度 较 高 。 


4 结论 


本 文 基 于 Landsat 8 OLI 影 像 的 时 序 特 征 ,构建 
信息 提取 方法 对 地 表 履 被 信息 进行 提取 。 首 先 根 
据 地 物 的 纹理 特征 和 旱 康 地 区 植被 的 物候 信息 选 
择 典 型 地 物 样 点 ,分析 地 表 罗 被 时 序 波谱 NDVI Al 
EVI 特 征 ,建立 典型 地 物 时 序 特 征 向 量 。 根 据 特 征 
向 量 匹配 和 最 近邻 原则 自动 确定 像 元 的 地 表 团 被 
类 别 ,完成 地 表 履 被 分 类 图 。 最 后 ,通过 新 疆 生 态 
与 地 理 研 究 所 野外 实地 考察 和 高 分 辨 率 卫 星 影像 
对 随机 选取 的 2500 个 样 点 数据 建立 混 清 矩阵 对 其 
进行 精度 评价 ,其 总 体 分 类 精度 为 97.2% ,Kappa 系 
数 达 0.9655 ,信息 捕获 较为 准确 ,分 类 效果 良好 。 主 
要 结论 如 下 : 

(1) Landsat 时 序 光 谱 特 征 可 较 好 的 反映 新 疆 
旱 康 地 区 农作物 的 季 相 节律 性 和 物候 变化 规律 万 


93.75 95.67 98.26 94.90 100 


Kappa 系 数 =0.9655 


其 捕获 农作物 在 不 同 生育 期 的 特征 ,为 地 物 识 别提 
供 多 维 资料 , 且 对 复 播 作 物 的 识别 ,潜在 提供 了 提 
高 识别 精度 的 途径 ,降低 “异物 同 谱 " 出 现 的 概率 ， 
增加 地 物 间 的 可 分 离 性 。 

(2) 新 疆 旱 康 地 区 是 典型 的 绿洲 农业 ,本 文通 
过 遥感 影像 时 序 特征 匹配 的 自动 识别 方法 ,为 及 时 
更 新 农作物 种 植 结构 提供 基础 积累 ,对 探索 农田 水 
循环 及 节 水 灌溉 措施 具有 显示 意义 。 

(3) 本 文通 过 时 序 特征 匹配 原则 对 有 限 范围 典 
型 地 表 禾 被 识别 进行 探讨 ,在 时 序 禹 感 影像 因 云 或 
气 溶胶 而 缺失 时 , 仍 需 继续 加 强 方法 研究 ,推演 更 
多 典型 植被 种 类 的 时 序 特征 ,为 在 更 广泛 的 区 域 推 
广 方法 打 基 础 。 
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Abstract: Surface land cover is a special landscape feature modeled by natural and anthropogenic activities, 
significantly influencing global or local climate, hydrothermal circulation, material transport, and diversity of 
terrestrial ecosystem communities. Although numerous global or regional land cover datasets have been published, 
challenges exist in the retrieval of detailed surface cover information, especially on the automatic identification of 
crop types. Furthermore, the classification accuracy of land-use and land cover maps needs improvement, specifically 
on crops in the agricultural management community. Using temporal Landsat 8 OLI multispectral images, we analyzed 
the spectral features of several dominant land cover types from temporal remote sensing images, with a specific 
emphasis on the unique phenological rhythm of individual plants and investigated the automatic retrieval of land 
cover types. We organized temporal Landsat 8 OLI surface reflectance data covering Fukang City of Xinjiang, 
China in the 2016 vegetation growth season. According to the Fukang statistical yearbook, prior knowledge of crop 
phenology, and the visible texture transition of vegetation in the growing season, we systematically selected the 
training samples of typical surface covers. We investigated the spectrum, the normalized difference vegetation in- 
dex, and the enhanced vegetation index temporal curves of 10 typical ground objects, such as agriculture fields of 
wheat, corn, and cotton, to extract the feature vectors modulated by the growth phase of each crop type. The Eu- 
clidean distance was calculated to measure the similarity between unknown surface targets and apriori templates 
built-up using the above feature vectors and assigned the land cover type according to the K-nearest neighbor 
tule. As a case study, we applied this approach to the land cover identification of Fukang using the time series of 
Landsat spectral data in 2018. We validated the retrieval accuracy using 2500 ground samples randomly collected 
from the land cover map of the investigated area, followed by the manual interpretation of the land cover types us- 
ing the high-resolution images of 2018 and our field surveys. The results correlate well with the overall accuracy 
of 97.2% and the Kappa coefficient of 0.9655. Specifically, combined with the field investigation results, the meth- 
od provides multi-dimensional data for ground object identification and potentially outlines the growth and develop- 
ment progress of different crops, enriches classification information, effectively reduces the phenomenon of for- 
eign matters with the same spectrum in a single time phase, and has potential advantages in improving the identifica- 
tion accuracy of re-plant fields. The time series spectral information captured by this research well reflects the sea- 
sonal rhythm and phenological change law of crops in the Fukang region of Xinjiang and provides a novel techni- 
cal method for the needs of agricultural monitoring. 

Key words: time series remote sensing data; phenology of vegetation; spectral characteristics; normalized difference 
vegetation index; enhanced vegetation index; surface cover mapping 


